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FOREWORD 


This report was prepared by the Pratt & Whitney Aircraft Division of 
United Aircraft Corporation, West Palm Beach, Florida, to describe the 
aerodynamic and mechanical design work conducted under Contract NAS3-11158, 
Task III, Single-Stage Experimental Evaluation of Tandem-Airfoil Rotor and 
Stator Blading for Compressors. Mr. Everett E. Bailey, NASA-Lewis Research 
Center, Fluid System Components Division, was Project Manager. 

The requirements of NASA Policy Directive NPD 2220.4 (September 14, 
1970) regarding the use of SI Units have been waived in accordance with the 
provisions of paragraph 5d of that Directive by the Director of Lewis Rearch 
Center. 
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summary 


Two compressor stages were designed for a comparative experimental 
evaluation in a 0.8 hub/tip ratio single-stage compressor. One stage was 
comprised of a conventional rotor and stator, and the other stage was comprised 
of a tandem rotor (two airfoils in series) and a tandem stator. The conventional 
and tandem rotors were both designed to produce a pressure ratio of 1.28 at a 
rotor tip velocity of 757 ft/sec. The stage pressure ratio at design flow and 
rotor speed was 1.26 for both stages. Predicted rotor and stage adiabatic 
efficiencies at design flow and rotor speed were 90% and 85%, respectively. 

A preliminary design study was performed to examine the effects of work 
input level and distribution on the rotor and stator diffusion factors, stage exit 
velocity profile, and overall efficiency. Based on this study, a radially constant 
work input distribution and an average rotor pressure ratio of 1.28 were 
selected for the rotor designs. Velocity diagrams and blade leading and trailing 
edge angles selected for the conventional rotor and stator blading were used in 
the design of the tandem blading. The effects of axial velocity ratio and secondary 
flow on turning were included in the selection of blade leading and trailing edge 

angles. 

Stress analyses were performed for the selected blading. These included 
analysis of blade attachment and disk stresses, vibratory stresses, and flutter. 
Materials that would provide adequate stress margins were selected for blade 
fabrication. 
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INTRODUCTION 


The effectiveness of tandem airfoils as a means for increasing the loading 
limit and stable operating range of highly loaded compressor blade rows was 
investigated for the National Aeronautics and Space Administration at the 
Florida Research and Development Center of Pratt & Whitney Aircraft under 
Task I of Contract NAS3-11158 (References 1 through 3). During this program, 
tandem rotors demonstrated higher pressure rise and efficiency than a con- 
ventional single airfoil rotor with identical inlet and exit airfoil angles. The 
performance of the conventional stage was controlled to a large extent by 
three-dimensional flow effects associated with high losses near the walls. The 
three-dimensional flows resulted even though the blading was designed with 
increased work input near the walls to compensate for the high losses in these 
regions and, thereby, maintain a constant radial pressure distribution. In 
principle, tandem rotors improved the performance by distributing the overall 
aerodynamic loading between the airfoils in tandem and effectively reduced the 
three-dimensional flows near the walls that were present with the highly loaded 
conventional rotor. 

Based on these results, a single-stage compressor investigation was 
initiated to evaluate the potential of tandem blading for improving the performance 
of a conventional stage designed with a selected radial work distribution that had 
lower work input at the walls. The radial work gradient was selected on the basis 
of rotor and stator diffusion factors in the end wall regions, stage exit flow 
profiles, and stage efficiency. The lower work input near the walls should 
reduce the three-dimensional flows and high wall losses that are characteristic 
of highly loaded blade rows and provide a comparison between a tandem stage 
and a conventional stage that is not characterized by a highly three-dimensional 
flow and associated poor performance. The aerodynamic and mechanical design 
of a conventional and a tandem airfoil stage is the subject of this report. The 
conventional single airfoil rotor and stator were designated Rotor D and Stator D. 
The tandem-blade rotor and stator were designated Rotor E and Stator E. 
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AERODYNAMIC DESIGN 


Design Guidelines 


The selection of the design velocity diagrams was accomplished within 
the range of the design guidelines given in table I. 

Table I. Design Guidelines 


Rotor Tip Diameter 

30 in. (Minimum) 

Hub/Tip Ratio 

0.7 to 0.8 

Rotor Tip Speed 

800 ft/sec (Maximum) 

Rotor Tip Diffusion Factor 

Less than 0.55 

Rotor Tip Solidity 

1.4 to 1.5 

Stator Hub Diffusion Factor 

Less than 0.60 

Stator Hub Solidity 

1.5 or greater 


In addition to the guidelines specified in table I, the following criteria 
were specified for the design: 

1. No inlet guide vanes (axial inlet flow) 

2. Axial discharge flow at the stator exit 

3. Common flowpath for both stages (same as used for 
Reference 1 program) 

4. Double-circular-arc blade sections. 

To ensure a valid comparison between the conventional Stage D and the 
tandem blade Stage E, the velocity diagrams selected for Rotor and Stator D 
were used to design the tandem blading. 
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SELECTION OF DESIGN PRESSURE RATIO AND ROTOR 
WORK DISTRIBUTION 


A preliminary design study of the effect of the radial distribution of work 
input on rotor diffusion factor, stator diffusion factor, stage exit flow profile, 
and efficiency was conducted to select the design rotor work distribution and 
pressure ratio for Stages D and E. Velocity diagrams were calculated for rotor 
pressure ratios of 1.28 and 1.32 for the following three radial distributions of 
work input: (1) constant spanwise, (2) that required to produce a constant rotor 
exit total pressure, and (3) parabolic (i.e., reduced work input near the walls). 
Velocity diagrams were also calculated for constant spanwise rotor work input 
at pressure ratios of 1.20 and 1.24. 

The initial phase of the design study involved an updating of the loss 
parameter vs diffusion factor correlation presented in Reference 1. To better 
define the radial loss profile, the loss correlation was expanded from three 
span locations (10, 50, and 90%) to five span locations (10, 30, 50, 70, and 90%). 
The 30 and 70% span data were obtained from the same data references that 
were used in the three-span loss correlation (Reference 1). The five-span loss 
correlation was then updated by adding data from NASA-sponsored programs 
(References 3 through 10) and unpublished loss data from in-house single-stage 
compressor programs performed at FRDC. The design loss curves selected 
to represent the data at each percent span are shown in figures 1 through 10. 

The two-dimensional cascade data from figure 149 of Reference 11 and the range 
of compressor data shown in figure 192 of Reference 11 are included for com- 
parison with the selected loss curves at 10, 50 and 90% span from the tip. 

The velocity diagrams were calculated by means of an iteration procedure 
using an axisymmetric flow field calculation and the loss correlations shown 
in figures 1 through 10. The calculation procedure solved the continuity, energy, 
and radial equilibrium equations, which included the effects of streamline 
curvature and radial gradients of enthalpy and entropy. The Reference 1 flowpath 
dimensions (figure 11), a rotor tip speed of 757 ft/sec, and an equivalent flow of 
110 lb/sec were maintained for the velocity diagram calculations. These values 
of rotor tip speed and flow are consistent with the design values used in the 
Reference 1 program and result in a rotor tip inlet Mach number of approximately 
0.8 and a specific flow of 33 lb/sec-ft 2 ; these values are generally representative 
of current design practice for compressor middle stages. 

The rotor inlet total pressure distribution from the data of the Reference 1 
program was used for the velocity diagram calculations. The resulting inlet 
velocity distribution is shown in figure 12. This velocity distribution represents 
an equivalent blockage of 3.0% of the total annulus area, and no additional 
blockage allowance was used at the rotor inlet. Figures 13 and 14 illustrate 
the differences in rotor inlet relative air angle and rotor and stator diffusion 
factors, respectively, which are associated with the measured rotor inlet total 
pressure (i.e., velocity) and a constant rotor inlet total pressure for an average 
rotor pressure ratio of 1.28 and constant rotor work. If the higher relative 
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air angles near the walls for the measured total pressure (i.e., velocity) profile 
are ignored during the metal geometry selection, the blade sections near the 
walls will be operating at high incidence angles. As discussed in Reference 12, 
these high incidence angles can stall the blade sections near the walls and affect 
the blade row performance over a large portion of the blade span. As shown in 
figure 14, ignoring the inlet total pressure gradients near the walls also results 
in substantial errors in the estimated values of rotor and stator diffusion factors 
near the walls. A blockage allowance of 5% of the local annulus area was 
assumed at both the rotor exit and stator exit to account for boundary layer 
growth on the flowpath walls at these locations. 

The resulting radial distributions of rotor diffusion factor, stator diffusion 
factor, and stator exit axial velocity are shown in figures 15 and 16 for the three 
radial distributions of work input investigated for rotor pressure ratios of 1.28 
and 1.32. Rotor diffusion factor, stator diffusion factor, and stage exit axial 
velocity are shown in figure 17 for the three pressure ratios (i.e., 1.20, 1.24, 
and 1.28) investigated for constant spanwise rotor work input. The stator tip 
region diffusion factor and exit axial velocity are not shown in figure 16 for the 
parabolic rotor work condition at a rotor pressure ratio of 1.32 because as the 
tip region loading was increased, the total pressure losses increased at a faster 
rate, and the iteration between the axisymmetric flow field calculation and the 
loss correlation would not converge. The overall performance in terms of 
pressure ratio and efficiency for the eight combinations of pressure ratio and 
radial distribution of rotor work input are summarized in table II. 


Table II. 

Summary of Predicted Overall Performance 


Work Distribution 


Rotor 

Stage 


Pressure 

Ratio 

Adiabatic 

Efficiency 

Pressure 

Ratio 

Adiabatic 

Efficiency 

Constant Work 

1.204 

91.6 

1.196 

88.3 

Constant Work 

1.246 

90.5 

1.235 

86.6 

Constant Work 

1.282 

89.9 

1.265 

84.8 

Constant Pressure 

1.295 

88.1 

1.272 

81.1 

Parabolic Work 

1.283 

91.9 

1.271 

87.6 

Constant Work 

1.319 

88.6 

1.294 

82.5 

Constant Pressure 

1.318 

87.3 

1.290 

80.2 

Parabolic Work 

1.323 

88.8 

* 

* 


^Values not available because the axisymmetric flow field calculation and 
loss correlation would not converge in the tip region as discussed above. 
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The radial diffusion factor distributions shown in figures 15 through 17 
indicate that too much rotor work input near the walls causes high rotor diffusion 
factors, whereas, too little work input results in high stator diffusion factors 
and a very nonuniform stage exit axial velocity (and flow). Although the para- 
bolic rotor work profile design resulted in the highest predicted stage efficiency, 
the stator end wall diffusion factors, shown in figures 15 and 16, were con- 
siderably higher than would be considered acceptable for design. Based on these 
results and the rotor and stator diffusion factor levels specified in table I, an 
average rotor pressure ratio of 1. 28 with a constant radial distribution of work 
input was selected for the Stages D and E design. These design conditions re- 
sult in a rotor diffusion factor of 0.50 at 10% span from the tip and a stator 
diffusion factor of 0.45 at 90% span, both of which are well within the respec- 
tive maximum allowable values specified in the design guidelines shown in 
table I. The radial distributions of stage inlet and exit velocity for the selected 
rotor pressure ratio of 1. 28 with constant rotor work input are shown in fig- 
ure 18. As shown in figure 18, the radial distribution of the stage exit velocity 
is slightly more nonuniform than the inlet velocity distribution. However, a 
secondary objective of this program is to evaluate the impact that a stage de- 
signed with constant work input would have on other stages. 

The velocity diagram calculation results selected for the Stage D and E 
design are shown in tables III and IV for the rotors and stators, respectively, 
along streamlines that pass through 5, 10, 15, 30, 50, 70, 85, 90, and 95% span 
at the rotor exit. The diffusion factor, loss coefficient, and exit total pressure 
distributions are also presented in figures 19 through 21. The predicted average 
rotor pressure ratio and adiabatic efficiency are 1.28 and 89.8%, respectively, 
at a rotor tip speed of 757 ft/sec. The predicted average pressure ratio and 
efficiency for the stage at design rotor speed are 1.26 and 84.8%, respectively. 

METAL GEOMETRY SELECTION 
Stage D 


Simulated double-circular-arc airfoils (i.e., the mean camber line and 
the suction and pressure surface lines of each blade element are lines with a 
constant rate of angle change with path distance on a specified conical surface), 
having constant chord length were selected for the rotor and stator. The 
thickness-to-chord ratio distributions, chord lengths, number of blades, and 
the number of vanes used in the design of the Reference 1 blading were used 
in the design of the Stage D blading. (See pages 12 and 14.) 

A study performed by Pratt & Whitney Aircraft has revealed better 
agreement between predicted and measured rotor and stator exit air angles 
when the cascade turning is modified to include the effects of axial velocity 
ratio and secondary flow. The predicted values both with and without the 
corrections for axial velocity ratio and secondary flow for the Reference 3 
blading are compared with the measured values in figures 22 and 23. Therefore, 
this technique was used in combination with equations 286 and 288 presented 
in Reference 11 to select the Stage D camber and incidence angles. 
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Table in. Rotor Vector Diagram Calculation Results 
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Table IV. Stator Vector Diagram Calculation Results 



CD D- CD CM CO O CD ^ CO 

onncot-coocoH^ 

• •••••••• 

OOOHOJ^^iOLO 

I I I I I 1 1 I I 


ooooooooo 

ooooooooo 


ooooooooo 

ooooooooo 


r-iCMCOOCMOOCMCM*^ 
• •••••••• 

lO(M^C 0 ^l>CDi >00 

Ot-H^lO^OOHD 

CO^lOOlOlO^^(M 


• ••••*••• 

lOCM^CO” 3 <t>Ot >00 

CO-^LOIDIDIO^^M 


COhOOHHHDCO't 

^Hooio^mioio 




t-aooooioocoOT— i 
• •••••»•• 
t-Nt-COON^C 5 ^ 
HHOCiCOCOin^^ 
Tt<T^Tj<COCOCOCOCOCO 


OOH<Mt>OOOCMH 
• •••••••• 

COOHcOWOHb-CO 
OOlOONCON 1>00 
CO^IDLOLOID^^W 


OOOCDOOOlOCMtHCM 

CDCDLOOLOCMlOOOO 

0 H^[OlO^O 5 COW 

LOCDCOCOCDCOLOLOt^ 


ooooooooo 


loolooooloolo 

DDCOOinc^HH 


ooooooooo 


LOOLOOOOLOOLO 
CD O 00 O- LO CO H H 




ajtxNoo^ww^a) 

HHt-^COOlOCOCO 


t>00 00 00 00 00 00!>CD 


CO MN ^ o CO Ifl lO W 
H(mohh 0 (oi>co 

NOClt-OCOHiiJN 

ooooooooo 

• •#•••••• 

ooooooooo 


WCOM^^t'LDCD^ 

t-Ot-HOCOOOOOCO 

OQOt^CDLOLOOLOCD 

OOOOOOOtHtH 

ooooooooo 


ONCOOO^OOOOO 

^CDNOOCOtDHOOCO 

• •••••••• 

ooooooooo 


^MHMNONOOOO 

CMCDLOCDCOOCD^LO 

O^M»OOMNl>C3 

LOLOLOlOLOLOlO^fCO 

• •••••••• 

ooooooooo 


ooooooooo 
• •••••••• 

LOOLOOOOLOOLO 

OOOOt>tOCOiHH 


ooooooooo 

LOOlOOOOLOOLO 

OftOOt-lOCOHH 



8 



For Rotor D, the turning that combines with the axial velocity ratio and 
secondary flow corrections to produce the required work input was calculated 
using the following procedure: 

1. A spanwise distribution of the required axisymmetric two- 
dimensional rotor exit relative air angle (/3j- e )2D was assume ^. 

The initial spanwise distribution was selected to be identical 
to the design rotor exit air angle (/3te) distribution shown in 
table III, and subsequent selections were based on the magnitude 
of the corrections calculated in step 4, below. 

2. The spanwise distribution of the axial velocity ratio (AVR) was 
calculated using the axisymmetric flow field calculation routine 
and the loss correlations shown in figures 1 through 10. 

3. The secondary velocities due to the blade row exit streamwise 
vorticity and the trailing edge filament vorticity were calculated. 
The streamwise vorticity, which includes the distributed 
passage and trailing shed vorticities, was calculated using 

the equation for the rate of change of the vorticity given in 
Reference 13: 

i.e. , V’ . V « ) * v ] • [< v ' • v > V ] + |r 2 [ n • (V X £ ) 


4. 


This equation was integrated to determine the exit streamwise 
vorticity by making the simplifying assumption that the tangential 
velocity varied linearly through the blade row. The trailing 
edge vorticity was calculated using the method of Reference 14. 


The assumed values of (£te)2D were corrected for the axial 
velocity ratio and secondary flow effects to obtain values of 

Pte * i • e • » 

/3 te = arctan — tan (j3 tea ) + V 0/ V zt e 

L Vz te J 

/3 t e a = arctan ^(AVR) a tan Ote^o] 


Where: 


V Zte = exit axial velocity (V zte + v z ) 


z te 


= axisymmetric exit velocity (i.e., flow 
a redistribution due to total pressure losses) 


(AVR) a = axisymmetric axial velocity ratio 


V0 = gap-averaged tangential secondary velocity 


v z = axial component of secondary velocity 
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5. The calculated values of p te were compared to the design 
velocity diagram angles,, If the calculated values of /?te did 
not agree with the design velocity diagram angles, a new 
spanwise distribution of the axisymmetric two-dimensional 
rotor exit air angle (Pte)2D was assumed and the procedure 
repeated. 

6„ The final values of (j3te)2D were then used, along with the 

known inlet conditions, to calculate the design two-dimensional 
turning for use in equations 286 and 288 of Reference 11 to 
calculate the incidence and camber angles. The three- 
dimensional connections given in Reference 11 were omitted 
since the selected turning including the effects of axial velocity 
ratio and secondary flow. Radial distributions of the final 
values for the following items are summarized in table V : 


1. 

Axial secondary velocity at the rotor exit (v z ) 

2. 

Gap-averaged tangential secondary velocity at the rotor 
exit (v 0 ) 

3. 

Axisymmetric two-dimensional rotor exit relative air 

angie (/3 te ) 2D 

4. 

The correction to the rotor exit relative air angle due 
to secondary flow (A/3 te )gp 

5. 

The correction to the rotor exit relative air angle due 
to the axial velocity ratio (A/^avr 

6. 

Rotor exit relative air angle corrected for secondary flow 
and axial velocity ratio (/?j- e ). 

Table V. 

Effect of Axial Velocity Ratio and Secondary Flow 
Corrections on Rotor Exit Air Angle 



Trailing Edge 








Percent Span 
From Tip 

v z 

v 9 

(0te)2D 


(A/3 te ) A vR 

Pte 

Hub 

95.0 

9.95 

0.95 

22.40 

-0.49 

4.64 

26.55 


90.0 

9.45 

1.35 

22.35 

-0.29 

1.24 

23.30 


85.0 

8.20 

2.40 

22.75 

-0.12 

0.37 

23.00 


70.0 

5.60 

2.50 

26.25 

-0.02 

-0.68 

25.55 


50.0 

2.55 

0.55 

30.55 

-0.08 

-1.27 

29.20 


30.0 

-1.10 

-3.40 

34.65 

-0.19 

-0.86 

33.60 


15.0 

-7.55 

-7.10 

38.45 

-0.09 

0.69 

39.05 


10.0 

-15.30 

-3.70 

40.90 

0.79 

1.31 

43.00 

Tip 

5.0 

-29.10 

10.50 

45.10 

3.25 

2.35 

50.70 
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The final calculated radial distributions of the axial velocity ratio and the rotor 
exit relative air angle (/3te) are compared with the design velocity diagram 
values in figure 24. The simulated double-circular-arc airfoil sections selected 
for the rotor were positioned on conic surfaces that approximate the design 
streamlines of revolution, and the resulting Rotor D metal geometry on these 
conic surfaces is summarized in table VI. The radial distributions of camber, 
incidence, and deviation angles are shown in figures 25 and 26. The design 
velocity diagram turning given in table III and the camber and incidence angles 
given in table VI were used to calculate the rotor deviation angle (i.e., 

5° = </> - A/3 + i m ). 


The following procedure was used to select the stator metal geometry: 

1. The secondary velocities due to the stator exit streamwise 
vorticity were calculated. 

2. The stator exit air angles (/3 £ e ) were modified to include the 
effect of secondary flow on turning to obtain the two-dimensional 
exit air angle that would theoretically result in an axial exit 
flow (i.e. , d te = arctan V(9 /V zt e ) * 

3. The radial distributions of the gap-averaged tangential 
secondary velocity (vg) and two-dimensional exit air angle 
(0te)2D are summarized in table VEL The camber angles 
required to turn the flow from the design velocity diagram 
inl et angle to the stator exit air angle calculated in item 2 
were calculated using equation 288 of Reference 11, except 
that the three-dimensional corrections for incidence and 
deviation angles were omitted. 

4. The camber angles required to turn the flow from the design 
velocity diagram inlet angle to the design exit angle (i.e., 0 deg) 
were calculated using the same method used in item 3. 

5. The actual design camber angles were obtained by averaging 
the camber angles calculated in items 3 and 4. (See figure 27.) 
This method was arbitrarily selected, since the predicted 
secondary flow resulted in more overturning near the walls 
than has been observed from the results of the test programs 
reported in References 3, 4, 5, 9, and 10. 

6. The incidence angles were calculated using equation 286 of 
Reference 11, except that the three-dimensional correction 
was omitted. 

7. The deviation angles were calculated from the design velocity 
diagram turning (table IV) and the camber and incidence angles 
calculated in items 5 and 6, respectively (i.e., 5° =<t> - A/3 + i m ). 

The resulting Stator D metal geometry along design streamlines is 
summarized in table VIII, and the radial distributions of camber, incidence, 
and deviation angles are shown in figures 28 and 29. 


11 



Airfoil: Simulated Double- No. of Blades: 70 Chord Length: 2.57 in. 
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Table VII. Stator Gap-Averaged Tangential Secondary Velocity 
and Two-Dimensional Exit Air Angle 



Trailing Edge 
Percent Span 
From Tip 

(0te)2D 



Hub 

95.0 

11.10 

-48.5 



90.0 

1.00 

-11.5 



85.0 

-0.75 

7.1 



70.0 

-0.95 

8.9 



50.0 

-0.80 

7.7 



30.0 

-1.35 

13.0 



15.0 

-0.35 

3.0 



10.0 

2.00 

-25.2 


Tip 

5.0 

22.70 

-140.5 


Rotor D and Stator D airfoil sections were designed on conic surfaces 
which approximate streamlines of revolution. For manufacturing purposes, the 
airfoil coordinates for each blade row were given on planes tangent to cylindri- 
cal surfaces normal to a radial line termed the stacking line. The airfoil sec- 
tions were positioned so that the stacking line passed through the center of 
gravity of the sections. A computer program provided a smooth fit of the air- 
foil properties and produced a set of coordinates for manufacturing purposes. 
Coordinates for the redefined sections, which were used for manufacturing 
purposes, are given in Appendix B. 


Stage E 

Tandem Rotor E and tandem Stator E were designed for approximately 
an equal distribution of loading between the front and rear airfoils; loading is 
defined as the tangential lift produced by the airfoil. To ensure interchange- 
ability with Stage D, the radial distribution of overall axial chord for the tandem 
rotor and stator blading were maintained equal to the values selected for the 
Stage D blading. To minimize the number of variables to be investigated in the 
selection of metal geometry, the maximum thickness-to-chord ratio for each of 
the airfoils in tandem was also maintained equal to the corresponding value 
selected for the Stage D blading. Simulated double-circular-arc sections were 
selected for both airfoils of the tandem blading. These airfoil sections were 
selected to closely approximate the type of blade sections used in an in-house 
study conducted at NASA-Lewis Research Center (Reference 15). The individual 
chords for the tandem blades were arbitrarily set equal, and the following pro- 
cedure was used to select the metal geometry. (This procedure was used in the 
Reference 1 design): 

1. Initial values of camber angle were assumed for the front and 
rear airfoils. The assumed camber angles were selected so 
that the passage between the blades would be slightly convergent 
(inlet-to-exit area ratio greater than one) to avoid undue velocity 
peaks or decelerations in the passage between the airfoils. 
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2. The individual airfoils were positioned according to the follow- 
ing criteria: 

a. The leading edge metal angle of the front airfoil and the 
trailing edge metal angle for the rear airfoil were main- 
tained equal to the leading and trailing edge metal angles, 
respectively, selected for Stage D. 

b. The passage width between the blades was maintained at 
approximately 10% of the individual airfoil chord. This 
selection was based on the results of the NASA in-house 
analytical study of tandem blading described in Refer- 
ence 15. 

c. Zero axial overlap of the front and rear airfoils was main- 
tained for ease of fabrication; however, this selection was 
consistent with the cascade results presented in Refer- 
ence 16 and the rotor results in Reference 17, and yielded 
a blade passage area ratio in the same range indicated 

as favorable in the NASA studies (Reference 15). 

3. Blade surface pressure and velocity distributions were calculated 
for two-dimensional, incompressible, inviscid potential flow. 

For these calculations, the airfoil sections on the conic surfaces 
that approximate design streamlines of revolution were assumed 
to be double-circular-arc sections on planes tangent to the 
conic and rotated to planes tangent to a cylindrical surface. 

The axis of rotation was defined by the intersection of two 
planes: one plane tangent to the cylindrical surface and the 
second plane normal to the compressor centerline. This second 
plane was located midway between the blade row leading and 
trailing edges . The potential flow solution involved a computer 
program that calculated the velocity field of an infinite cascade 
as governed by Laplace's equation: 

2 

V <t> = 0 

where 4> is the velocity potential. Solutions were obtained for 
a zero angle of attack, 90-deg angle of attack, and circula- 
tory flow for each airfoil, and the results superimposed in 
such a way that the design angle of attack was obtained and the 
Kutta condition satisfied. The method of solution, described 
in Reference 18, uses a distribution of sources on the airfoil 
surface and solves a set of linear algebraic equations for the 
source distribution that forces the total velocity normal to the 
airfoil surface to be zero. The total velocity is the sum of two 
velocities: the onset velocity, defined as the velocity field in 
which the body is immersed, and the disturbance velocity, 
defined as the velocity field caused by the source distribution. 
The blade surface pressure and velocity distributions thus 
obtained were corrected for compressibility by means of the 
Karman-Tsien equation. 
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4. The difference in the pressure surface and suction surface 
pressures were integrated to determine the loading (tangen- 
tial lift) for each of the two airfoils in tandem. 

5. The maximum suction surface-to-exit velocity ratio for each 
airfoil was calculated to ensure that it was not greater than 
2.0. According to the results presented in Reference 19, a 
velocity ratio in excess of 2. 0 may lead to a possible rapid 
increase in the airfoil boundary layer momentum thickness 
with a corresponding increase in loss. The maximum velocity 
on the suction surface was increased to reflect streamtube 
convergence and radius change through the blade passage 
that were not accounted for in the two-dimensional analysis. 

The amount of increase was the local difference between 
linear distributions of velocity from the inlet values to the 
exit values obtained from (1) the velocity diagram calculations 
and (2) the potential flow calculations, as illustrated by LC 

in figure 30. The velocity ratio of interest was the corrected 
suction surface velocity divided by the exit velocity obtained in 
the velocity diagram calculations. 

6. The individual camber angles were changed while holding the 
overall total camber angle constant and the procedure repeated 
until approximately an equal distribution of loading was obtained 
between the two airfoils in tandem and no sharp velocity peaks 
existed in the passage between the airfoils. 

The resulting spanwise camber distributions for Rotor E and Stator E 
are shown in figures 31 and 32, respectively, and the Stage E metal geometry 
is summarized in table EX. 

The static pressure coefficients and velocities at the radial locations 
used to define the spanwise camber angle distributions are shown in figures 33 
through 64 for the final tandem airfoil configurations. The maximum suction 
surface velocity corrected for streamtube convergence through the blade 
passage and the design inlet and exit velocities obtained from the velocity 
diagram calculations are shown on the appropriate figures. The loading splits, 
camber ratios, and the airfoil maximum suction surface-to-exit velocity ratios 
are given in table X. The blade passage overlaps, gaps, and convergences 
(figure 65) are also given in table X. 

The Stage E airfoil sections were designed on conic surfaces, and the 
airfoil coordinates were given on plane surfaces tangent to cylindrical surfaces 
for manufacturing purposes using the same procedure used for Stage D. The 
Rotor E airfoil sections were positioned so that the stacking line passes through 
the centers of gravity of the combined airfoil sections. The Stator E airfoil 
sections were positioned so that the stacking line passes through the trailing 
edges of the front airfoil sections. Coordinates for the tandfem airfoil sections 
which were used for manufacturing purposes are given in Appendix C. 
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Table X. Tandem Airfoil Design Summary 
Tandem Rotor E 
Work Split: 50% - 50% 
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* Information included in this table is defined on planes tangent to the conic surfaces which approximate design 
streamlines of revolution. 



MECHANICAL DESIGN 


Rotor Steady-State Stress Analysis 


Rotor D 

The stresses due to centrifugal loads and/or gas bending loads were cal- 
culated at 13 radial locations for a Rotor D airfoil fabricated from AMS 5616 
(stainless steel). The reduction of gas bending stress due to centrifugal load 
was considered, and the resulting net gas bending stress and centrifugal tensile 
stress were added to yield the total blade stress at each radial location. The 
results of this analysis are presented in figure 66. The maximum stress for 
Rotor D was 12, 000 psi in the leading and trailing edges of the hub. This cal- 
culated stress is well within the 0. 2% yield strength of 110,000 psi for AMS 5616. 

Rotor E 

Preliminary analysis of the front and rear airfoil natural frequencies for 
Rotor E indicated that a bridge connecting the two airfoils in tandem was re- 
quired to increase frequency and stiffness (thus reducing susceptibility to flutter, 
as will be discussed on page 21) and to ensure dimensional stability during 
operation. Rotor E, with a 0.060-in. thick interblade bridge at 30% span, was 
analyzed to determine airfoil and interblade bridge stress due to centrifugal and/ 
or gas bending loads. To minimize centrifugal force stress, titanium (AMS 4973) 
was selected as the blade material in preference to stainless steel (AMS 5616). 

An analytical model comprised of statically loaded elastic structures represented 
by slender prismatic beam members was used to determine tandem blade stress. 
Figure 67 presents a graphic description of the analytical model. The front and 
rear airfoils of Rotor E were each divided into 10 elements, and the interblade 
bridge was divided into 3 elements. The beam members or elements were repre- 
sented by their centroidal axis and analyzed as line elements. Centrifugal and 
aerodynamic loads were then applied to each element to yield reactions and dis- 
placements of the element or the element end, i.e., joint. These values were 
then used to calculate the total stress values shown in figure 68 for Rotor E. A 
maximum stress of 23,000 psi occurred in the leading edge of the front airfoil 
at the hub. The maximum interblade bridge stress was 17,800 psi. These cal- 
culated steady-state blade element stresses are well within the 0.2% yield 
strength of 104,000 psi for titanium (AMS 4973). 

Rotor Vibratory Analysis 


Rotor D 

Bending and torsional vibratory frequencies were calculated for Rotor D 
and the results presented in terms of frequency vs rotor speed in figure 69. At 
design equivalent rotor speed, the calculated bending and torsional vibratory 
frequencies were 440 Hz and 1370 Hz, respectively. Lines representing multiples 
of rotor frequency (E) are shown in figure 69 to permit identification of resonant 
operating conditions that might possibly be encountered during testing due to 
upstream bearing support struts, rotating stall zones, or upstream instrumenta- 
tion. As can be seen in figure 69, the first bending mode for Rotor D is close 
to the 6E line at 100 and 110% design equivalent rotor speed. This indicates a 
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potential resonance condition because there are six inlet struts. However, no 
resonance condition is expected because 6E excitation frequencies have not been 
observed during testing of similar blading in the same compressor. The lack 
of 6E excitation frequencies is attributed to the unequally spaced struts (t/c =0.12) 
being located three chord lengths upstream of the rotor, thus providing sufficient 
distance for the strut wakes to be substantially dissipated prior to reaching the 
rotor inlet. No torsional vibratory problem is expected because of the low energy 
associated with the high multiples of rotor speed, shown in figure 69, required 
to produce a resonant condition. The excitations due to rotating stall cells would 
be in the range of 35 to 70 Hz at design speed and are well below any resonance 
frequency. Because of the relatively high resonant frequencies shown in fig- 
ure 69 and the large vibratory stress margin available (as indicated by the 
Goodman diagram of figure 70, i.e., 56,000 psi based on the smooth fatigue 
strength and 19, 500 psi based on the notched fatigue strength), no Rotor D vibra- 
tory problem is anticipated. 

Rotor E 

Bending vibratory frequencies were calculated for the individual front and 
rear airfoils and the bridged airfoil configuration of tandem Rotor E and the 
results shown in terms of frequency vs rotor speed in figure 71. The bridged 
blade bending frequencies were determined by restraining the front and rear 
airfoils at the bridge location so that bending vibratory motion was permitted 
in one plane, i. e. , at one shroud angle. The shroud angle was then varied until 
both the front and rear airfoils achieved the same frequency, or theoretically 
vibrated together. The frequency so calculated for tandem Rotor E at design 
equivalent rotor speed was approximately 470 Hz. The Rotor E bending fre- 
quency approaches the 6E line at 110% design equivalent rotor speed, but, as 
with Rotor D, no resonance problem is anticipated for this rotor. The same tech- 
nique was used to calculate the bending frequencies of the Reference 2 tandem 
blading, and good agreement was obtained between the measured and calculated 
frequencies. 

Torsional vibratory frequencies calculated for the front and rear airfoils 
of Rotor E at design equivalent rotor speed were 740 and 770 Hz, respectively. 

A torsional frequency for the bridged blade configuration was not calculated. 
However, the bridged blade frequency is not expected to deviate substantially 
from the individual airfoil frequencies. When the frequencies for the Reference 2 
tandem blading were recalculated using the same technique used for Rotor E, 
the frequencies of the individual airfoils were within 11% of the value measured 
for the bridged blade configurations. No torsional resonance problem is ex- 
pected since, as shown in figure 71, the individual airfoil frequencies were well 
above the 6E excitation frequency. To illustrate the vibratory stress margin 
present in the design of Rotor E, a Goodman diagram for AMS 4973 is presented 
in figure 72. As shown in figure 72, allowable vibratory stress to failure for 
Rotor E was 43, 500 psi based on the smooth fatigue strength and 27,500 psi 
based on the notched fatigue strength. Neither configuration indicated a vibra- 
tory problem. 
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Rotor Flutter Analysis 


Rotor D 

Values of the reduced velocity and incidence parameters were calculated 
for Rotor D at the design operating conditions and for the estimated negative 
and positive incidence angle operating limits between 50 and 110% design equiva- 
lent rotor speed. These results were compared to correlated flutter data for 
the first bending and first torsional vibratory frequencies (figure 73). The re- 
duced velocity parameter is defined as: 

7TCU) 

and the incidence parameter is defined as: 

. _ Sn ~ * m ref 

^nr ~~ low-loss incidence range 

where V, c, im re f» anc * low-loss incidence range are the values for airfoil sec- 
tions located at 25% span from the tip. The low-loss incidence range and im re f 
were determined from an unpublished P&WA cascade data correlation. The 
bending and torsional mode flutter calculations were made at Mach numbers of 
0.4 and 0.6, respectively, so that the values obtained could be compared with 
the correlated data. No flutter problem is anticipated because, as indicated 
in figure 73, the operating range is in the safe operating zone. The operating 
region for the conventional airfoil of Reference 3, which did not encounter a 
flutter condition, is also shown in figure 73. Since the calculated operating re- 
gion for Rotor D compares favorably with the operating region for the Reference 3 
rotor, this further supports the prior conclusion that no bending or torsional 
flutter problems are anticipated for Rotor D. 

Rotor E 

Bending and torsional mode reduced velocity parameters of 2.80 and 1.76, 
respectively, were calculated for the front airfoil of Rotor E. Because of the 
high value of the reduced velocity parameter for the front airfoil, a 0. 060-in. 
thick interblade bridge was added to Rotor E at 30% span to increase the blade 
natural frequency. 

Bending and torsional mode reduced velocities and incidence parameters 
were calculated for the bridged configuration of Rotor E at the design operating 
point and for the estimated negative and positive incidence angle operating limits 
between 50 and 110% design equivalent rotor speed and the results compared to 
correlated flutter data in figure 74. The reduced velocities and incidence param- 
eters for the bridged blade were calculated using the overall chord dimensions 
and the front airfoil incidence angles and velocities. This was done because the 
front airfoil of the tandem configuration is subjected to incidence angle variations, 
while the incidence angle variations on the rear airfoil are expected to be small 
because of the small variations in exit air angle from the front airfoil. The 
overall chord was used because the bridged blades will move together in the 
immediate bridge region. The Rotor E reduced velocity parameters for the 
bending operating range were based on the calculated bridged blade frequency. 
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However, since the bridged blade torsional frequency was not available, the 
rear airfoil frequency was used to calculate the torsional reduced velocity param- 
eter. Because individual airfoil frequencies are expected to be less than the 
bridged blade frequency, any conclusion based on the individual airfoil should be 
conservative. Values of the reduced velocity and incidence parameter for tandem 
Rotor E at the design operating conditions and the estimated negative and posi- 
tive incidence limits are shown in relation to correlated flutter data in figure 74. 

No flutter problem is anticipated. The operating region for the Reference 2 tandem 
blade, which did not encounter flutter, is also shown on figure 74 for comparison 
with the calculated operating region of Rotor E. Since the operating region for 
Rotor E is similar to the operating region of the Reference 2 blading, this further 
supports the prior conclusion that no bending or torsional flutter problems should 
occur during testing of Rotor E. 


Rotor Attachment 

Blade spindle tensile, bending, shear, and bearing stresses were calculated 
considering the airfoil centrifugal forces and gas bending stresses due to aero- 
dynamic loading. The calculations were performed at a rotor speed of 6000 rpm, 
which is approximately 140% of design speed. The results of the stress calcula- 
tions for Rotors D and E are presented in table XI. The combined tensile and 
bending spindle stresses calculated for Rotor D and Rotor E were 67, 700 and 
30,200 psi, respectively. These calculated stresses are well within the 0.2% 
yield strengths of 114,000 and 104,000 psi of the stainless steel (AMS 5616) and 
titanium (AMS 4973) selected for Rotor D and Rotor E, respectively. Similarly, 
calculated shear and bearing stresses of 12,900 and 84,000 psi for Rotor D and 
5, 760 and 37,500 psi for Rotor E, respectively, did not exceed the specified 
material limitation shown in table XI (i. e. , allowable shear stress equals 55% 
of material ultimate tensile strength and allowable bearing stress equals 120% 
of material 0. 2% yield strength). Consequently, no blade attachment stress 
problems are anticipated. 


Rotor Disk and Carrier 


The average tangential stress for the AMS 6415 (low alloy steel) rotor disk 
and carrier was determined through the use of a computer disk analysis program 
and found to be well within design practice for AMS 6415 (0. 2% yield strength 
of 140, 000 psi). The results of this analysis are presented in table XII. 


Stator Steady-State Stress Analysis 


The gas bending stresses in the leading and trailing edge and on the con- 
cave surface at the point of maximum thickness were calculated for Stator D and 
the front and rear airfoils of Stator E. Calculations were made assuming (1) 
the vanes would be fabricated from AMS 5613, a stainless steel that has a 0.2% 
yield strength of 110, 000 psi and (2) the vanes were beams that would deflect as 
guided cantilevers about the tip. The guided cantilever about the tip condition was 
selected, even though the stator vanes are attached to the shrouds by trunnions 
at both the hub and tip, to provide a conservative estimate of the vane stresses. 


22 



Table XI. Blade Attachment Stress (6000 rpm) 
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Table XII. Disk and Carrier Stress (6000 rpm) 


Configuration 

Rotor D 

Rotor E 

Disk (AMS 6415) Stress, psi 

94,000 

42,000 

Carrier (AMS 6415) Stress, psi 

62,000 

27,700 

0. 2% Yield Strength of 



AMS 6415 at 100 °F, psi 

140,000 

140,000 


Some movement of the vane at the hub is possible. These stress values are 
shown in figure 75 for Stator D and in figure 76 for Stator E. The maximum 
bending stress of 2, 000 psi occurred in the leading edge tip of Stator D. A maxi- 
mum Stator E front airfoil bending stress of 7,300 psi (compressive) was cal- 
culated for the convex surface at the tip. The maximum rear airfoil stress was 
8,300 psi and occurred in the leading edge at the tip. None of the calculated 
stator stresses were prohibitive because of the high yield strength of the vane 
material and no stress problem is anticipated. 

Stator Vibratory Analysis 

Bending and torsional vibratory frequencies were calculated for Stator D 
and Stator E front and rear airfoils and the results presented in terms of fre- 
quency vs rotor speed in figures 77 and 78. The vibratory analysis was made 
assuming the stators to be beams with both ends fixed (fixed-fixed mode). This 
assumption was permissible because stator hub and tip trunnions are held in 
inner and outer diameter shrouds. Lines representing multiples of rotor passing 
frequency (70E) are shown in the figures to permit the identification of any ex- 
citation frequencies within the operating range. Stator D bending and torsional 
frequencies in the fixed-fixed mode were 4650 and 4050 Hz, respectively. 

Stator E front and rear airfoil bending frequencies were 1820 and 1870 Hz, re- 
spectively, while the torsional frequency of the front and rear airfoils were 
4200 and 4000 Hz, respectively. No vibratory fatigue problems are anticipated 
for either Stator D or Stator E. Because of the low steady-state stress present 
in these vanes, a large vibratory stress margin is available, as indicated by 
the Goodman diagram of figure 79 for Stator D and figure 80 for Stator E. As 
shown in figure 79, Stator D can withstand 56, 000 psi vibratory stress based on 
the smooth fatigue strength and 21, 000 psi based on the notched fatigue strength. 
Similarly shown in figure 80, both Stator E front and rear airfoils can withstand 
53,000 psi vibratory stress based on the smooth fatigue strength; and 20, 000 psi, 
based on the notched fatigue strength. 


Stator Flutter Analysis 


Stator D and Stator E front and rear airfoil torsional stall flutter charac- 
teristics were calculated and presented for comparison with correlated flutter 






data in figures 81 and 82. The flutter variables are a reduced velocity param- 
eter K, as defined in the rotor flutter analysis section, page 24 and an average 
row pressure ratio, defined as: 


Pie + Pfe 
2 

Pie 

The values of velocity and chord used to calculate the stator reduced 
velocity parameter are the values at 50% span. Torsional frequencies for the 
stators in a fixed-fixed configuration were used in the calculations. As shown 
in figures 81 and 82, no flutter problems are anticipated for Stator D or the 
front and rear airfoils of Stator E. 

Stator Attachment 


Stator assembly is achieved by tack welding the cylindrical trunnions at 
each end of the vane into the inner and outer diameter shrouds. The cross 
section of primary interest for stator stress evaluation is the junction of the 
airfoils and trunnion. For this cross sectional area of Stator D, the calculated 
bending stress due to aerodynamic loading was 6, 530 psi. For Stator E front 
and rear airfoils, this stress was 15, 700 psi and 14,600 psi. These trunnion- 
airfoil stress values are well within the 0. 2% yield strength of 110, 000 psi for 
the AMS 5613 stainless steel material selected for stator fabrication, and no 
stress problems are anticipated. 
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Figure 1. Rotor Loss Parameter vs Diffusion Factor, 10% Span From Tip DF 93388 












Figure 2. Rotor Loss Parameter vs Diffusion Factor, 30% Span From Tip DF 93389 















28 


Figure 3. Rotor Loss Parameter vs Diffusion Factor, 50% Span 
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Figure 4. Rotor Loss Parameter vs Diffusion Factor, 70% Span From Tip DF 93391 
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Figure 6. Stator Loss Parameter vs Diffusion Factor, 10% Span From Tip DF 93393 
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Figure 8. Stator Loss Parameter vs Diffusion Factor, 50% Span DF 93395 
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Figure 9. Stator Loss Parameter vs Diffusion Factor, 70% Span From Tip 
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Figure 10. Stator Loss Parameter vs Diffusion Factor, 90% Span From Tip DF 93397 
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Figure 12. Rotor Met Absolute Velocity Distribution DF 93398 
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Figure 14. Effect of Compressor Met Velocity on DF 93400 

Rotor and Stator Performance 
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Figure 15. Effect of Rotor Work Distribution on DF 93401 

Rotor Loading, Stator Loading, and Stage 
Exit Axial Velocity for Rotor Pressure 
Ratio of 1. 28 
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Figure 16. Effect of Rotor Work Distribution on Rotor DF 93402 
Loading, Stator Loading, and Stage Exit 
Axial Velocity for Rotor Pressure Ratio 
of 1.32 
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Figure 18. Stage Inlet and Exit Axial Velocity Distri- DF 93404 
butions for Stages D and E 
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DF 93405 
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Figure 19. Rotor and Stator Diffusion Factor Distri- 
butions for Stages D and E 







Figure 20. Rotor and Stator Loss Distributions for DF 93406 
Stages D and E 
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Figure 21. Rotor and Stator Exit Pressure Profiles DF 93407 
for Stages D and E 




Figure 22. Effect of Axial Velocity Ratio and Secon- 
dary Flow Corrections on Rotor Exit Air 


DF 93408 
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Figure 24. Predicted Values of Rotor Exit Air Angle 

and Axial Velocity Ratio for Stages D and E 


DF 93410 
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Figure 27. Comparison of Stator Camber Angle Dis- 
tributions 


DF 93413 
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Figure 28. Stators D and E Incidence and Deviation 
Angle Distributions 


DF 93414 
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Figure 30. Description of Technique Used to Modify DF 93416 
Two-Dimensional Potential Flow Solution 
for Streamtube Convergence Through the 
Blade Row 
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Figure 33. Rotor E Static Pressure Coefficient Dis- DF 93419 
tribution, 0% Span From Tip 
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Figure 34. Rotor E Blade Surface Velocities, DF 93420 

0% Span From Tip 
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Figure 35. Rotor E Static Pressure Coefficient Dis- 
tribution, 2. 5% Span From Tip 


DF 93421 
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Figure 36. Rotor E Blade Surface Velocities, 
2. 5% Span From Tip 


DF 93422 
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Figure 38. Rotor E Blade Surface Velocities, 
9.3% Span From Tip 


DF 93424 






Figure 39. Rotor E Static Pressure Coefficient Dis- DF 93425 
tribution, 12. 2% Span From Tip 
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Figure 43. Rotor E Static Pressure Coefficient 
Distribution, 66.5% Span From Tip 


DF 93429 
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Figure 44. Rotor E Blade Surface Velocities, 66.5% DF 93430 
Span From Tip 
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Figure 46. Rotor E Blade Surface Velocities, 
91.0% Span From Tip 


DF 93432 
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Figure 48. Stator E Vane Surface Velocities, 
2.5% Span From Tip 


DF 93434 
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Figure 49. Stator E Static Pressure Coefficient DF 93435 

Distribution, 5% Span From Tip 
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Figure 50. Stator E Vane Surface Velocities, DF 93436 

5% Span From Tip 
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Figure 51. Stator E Static Pressure Coefficient Dis- DF 93437 
tribution, 9.1% Span From Tip 
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Figure 52. Stator E Vane Surface Velocities, DF 93438 

9.1% Span From Tip 
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Figure 53. Stator E Static Pressure Coefficient 
Distribution, 20% Span From Tip 


DF 93439 
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DF 93441 


Figure 55. Stator E Static Pressure Coefficient Dis- 
tribution, 50% Span 
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Figure 57. Stator E Static Pressure Coefficient Dis- 
tribution, 80% Span From Tip 


DF 93443 
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Figure 59. Stator E Static Pressure Coefficient Dis- 
tribution, 91. 1% Span From Tip 


DF 93445 
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Figure 60. Stator E Vane Surface Velocities, DF 93446 

91.1% Span From Tip 
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Figure 61. Stator E Static Pressure Coefficient Dis- DF 93447 
tribution, 95% Span From Tip 
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Figure 62. Stator E Vane Surface Velocities, DF 93448 

95% Span From Tip 
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Figure 63. Stator E Static Pressure Coefficient Dis- 
tribution, 97.5% Span From Tip 


DF 93449 
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Figure 64. Stator E Blade Surface Velocities, DF 93450 

97.5% Span From. Tip 
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BLADE PASSAGE CONVERGENCE - F = H/G 
BLADE PASSAGE GAP RATIO - G/c 
BALDE PASSAGE OVERLAP RATIO - L/c 
c = INDIVIDUAL AIRFOIL CHORD 
c c = COMBINED OR OVERALL CHORD 
c cx = OVERALL AXIAL CHORD 

Figure 65. Tandem Airfoil Geometry, Simulated 
Double-Circular-Arc Airfoils 


FD 64418 
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Figure 66. Calculated Rotor D Stress Distribution DF 93451 
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Figure 67. Graphic Description of Tandem Rotor FD 62295A 

Analytical Model 


92 




93 


Figure 68. Calculated Rotor E Stress Distributions DF 93452 
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Figure 70. Rotor D Goodman Diagram DF 93454 
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Figure 71. Rotor E Resonance Diagram for First Bending and First Torsional Mode Vibration 
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Figure 74. Calculated Rotor E First Bending and First Torsional Mode Flutter Characteristics DF 93458 
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Figure 75. Calculated Stator D Stress Distribution DF 93459 


100 





101 







3 




102 


Figure 77. Stator D Resonance Diagram for First Bending and First Torsional Mode Vibration DF 93461 





103 














Figure 79. Stator D Goodman Diagram 


DF 93463 
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Figure 80. Stator E Goodman Diagram 


DF 93464 
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Figure 81. Calculated Stator D First Torsional Mode DF 93465 
Flutter Characteristics 
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Figure 82. Calculated Stator E First Torsional Mode Flutter Characteristics DF 93466 





APPENDIX A 


AVR 

a 0 

c 

C P 

d 

D 

E 

F 

f(im) 

G 

H 

i„_ 

m 

K 

Kip 

L 

LC 

M 

N 

P 

PR 

P 

t 

T 

U 

V 


DEFINITION OF SYMBOLS 
Axial Velocity Ratio (V zte /V Z j e ) 

Inlet relative stagnation velocity of sound, ft/sec 
Chord length, in. 

Static pressure coefficient 
Diameter, in. 

Diffusion factor 
Multiple of rotor frequency 

Tandem airfoil passage convergence (See figure 61. ) 

Incidence parameter 

Tandem airfoil passage gap (te), in. (See figure 61. ) 

Tandem airfoil passage gap (le), in. (See figure 61. ) 
Incidence angle, deg 
Reduced velocity parameter 
Notch factor 

Tandem airfoil passage overlap, in. (See figure 61. ) 

Local correction for suction surface velocity (See figure 26. ) 

Mach number 

Rotor speed, rpm 

Total pressure, psia 

Pressure ratio 

Static pressure, psia 

Blade maximum thickness, in. 

Total temperature, °R 
Rotor speed, ft/sec 
Velocity, ft/sec 



DEFINITION OF SYMBOLS (Continued) 


v 

W 

x 

y 

z 

a 

P 

A/3 

y 

y° 

8 

8 ° 

r 

(9 

K 

p 

a 

<t> 

4> 

a 

a> 

w 

cucosffte 


Secondary velocity, ft/sec 
Actual flowrate, lb m /sec 
Radial direction 
Tangential direction 
Axial direction 

Cone angle (angle of plane tangent to conic surface that 
approximates the design streamline of revolution), deg 

Air angle, degrees from axial direction 

Flow turning angle, deg 

Ratio of specific heats 

Blade-chord angle, degrees from axial direction 

Ratio of total pressure to NASA standard sea level pressure 
of 14.694 psia 

Deviation angle, deg 

Absolute vorticity 

Ratio of total temperature to NASA standard sea level tempera- 
ture of 518. 7 °R 

Blade metal angle, degrees from axial direction 
Absolute vorticity in direction of relative velocity 
Density, lbf - sec 2 /ft^ 

Solidity 

Velocity potential 
Blade camber angle, deg 
Angular velocity vector 
Frequency, Hertz 
Loss coefficient 
Loss parameter 
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DEFINITION OF SYMBOLS (Continued) 

Subscripts: 

a Axisymmetric calculation value 

AVR Axial velocity ratio 

b Bending 

c Combined or overall 

cx Overall axial 

fs Freestream value 

id Isentropic condition 

le Leading edge 

ref Minimum loss 

SF Secondary flow calculation value 

te Trailing edge 

t Torsional 

z Axial component 

6 Tangential component 

2D Two dimensional 

Superscripts: 

' Related to rotor blade 

Mass average value 

DEFINITION OF DESIGN VARIABLES 

Incidence Angle: 

Rotor: i m - (3 le - K le Stator: i m = /? le - * le 

Diffusion Factor: 

Rotor: 


D = 1 - 


V\ d, V, - d, V fll 
te te le A le 


(d le + d te )crV’ le 
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DEFINITION OF DESIGN VARIABLES (Continued) 


Stator: D = 1 


V te d te V 9te d le V »le 
V le"< d le +d te> ffV le 


Deviation Angle: 


Rotor: «‘= ?(e* K te 


Stator: 6 = ^ - * te 


Loss Coefficient: 


- - P ’ 


Rotor: u)' = 


where: 


te'id te 

P'i - P, 
le le 


< P ’te>id * P 'le 


1 + 


7-l/ U te 


2 \ 

U Oie / L 


7 


'le 


te/ 

7 


17 - 1 


P' is found from p/P' = 


Stator: u; = 


P f - P^ 
fs te 

P p - P, 
fs le 


Static Pressure Coefficient: 


1+ 2|1 m .2U-7 


c = 


P L “ p fs 


P 1/2 Pts V2 ts 
Pressure Ratio: 


Rotor: 

Equivalent Flow: 

Wi/F 


rotor te 


rotor le 


Stage : 


stator te 


rotor le 
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DEFINITION OF DESIGN VARIABLES (Continued) 


Equivalent Rotor Speed: 

n/VT 

Adiabatic Efficiency: 


Rotor: ±r 




lz± 

y 


- 1 


T. /518. 7-1 
te 


Stage : 


y - 1 

(PR) 7 - 1 

T\ /518. 7-1 
te 
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APPENDIX B 


STAGE D AIRFOIL COORDINATES 


NOMENCLATURE FOR TYPICAL 
VIEW OF AIRFOIL EDGE RADIUS 



GRAVITY) 

ROTOR D STACKING CONFIGURATION 



STATOR D STACKING CONFIGURATION 


Figure B-l. Stage D Airfoil Coordinates 


FD 64419 
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LEI< = C.0363 Tlk = 0.0038 LER = 0.0063 TER * 0.0038 LEK « 0.006/ TEK • 0.003y 

Chord = 2.5837 CHORD = 2.5828 CHORD « 2.5785 



Table B-l. Stage D Airfoil Coordinates (Continued) 
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t.E. CHICLE CENTER -C.880? -0.92C7 -0.9074 -0.8969 -0.9345 -0.8712 

T.E. CIRCLE CENTER 1.0545 0.7351 1.0817 0.6902 1.1076 0.6436 

LEK . 0.0065 TEN » 0.0063 LEK = 0.0070 TEK « 0.0069 LEn » 0.0U76 TtK ■ J.0076 

CHORD » 2.5597 CHORD * 2.5587 CHORD * 2.5578 



Table B-l. Stage D Airfoil Coordinates (Continued) 
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LE CENTER 1.137C 0.5890 1.1204 0.6236 1.0960 0.673T 

LEW . 0.0091 TEK • 0.0091 LER . 0.0092 TEW ■ 0,0091 LEW . 0.0093 TtW « 0.0092 

CHORD = 2.5590 Chord = 2,5614 Chord = 2.5u7d 
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Table C-l. Stage E Airfoil Coordinates 
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Table C-l. Stage E Airfoil Coordinates (Continued) 


! 


NO 

ao.m'irtir-ir-N 

00 ® iN-|® 'm > 

• •••:• 

o:o o;o o < 

1 I i I I 1 I 1 


• O — 4'®|4 ® 
> m — ®rvtn ® 
I ® >0 to (M O' 4 
\ rvj — * ©O O — 

no o o!o o o 

III 


z. ui ! 

'A 

io'®oc I— | 
to • d! I 


>i— in- <*>ioi® — ■ r- rui o ® 

v® ® r- ®!® O' o olO 1 o ® 

> mm to: to to to to; ®> — ® 

jn-® ®4mcg — oj O o — 

VO O o oio O O O' o o o 

I I I I I I I I 


!® -4: to ® O • 
| O' O'lmif*-®' 
> — f>- rgiinir- < 

B N N- ®itn • 
• • • • • 

0 o o o o « 

1 I I I I 


® 4 o 4 n i 

•to ® ® o o t 

• N- — « f«V4 <0 f 
® ® f- ® m - 


II® Olio 0<-0 O'; 
»< 4 l®;m| 0 'i(O rg 
i|r»llO!N“ O lO — 

■ r^t I rvi I • — © O' 

I • . 

) 0i0;0 O O O 

i;i ii i i : 


) co in O'iHWjn 
• in ® \0 to.o' 4 ■ 
j rg pg pg: 4 - o to 
»rg — OO — — 


O O OO O 0:0 o 0.0 OO O 

i . i III » ! i i . i ( 

. I 

'N-cg' — irgi — iOlN-:N“.® iO — O ® i 
ICtIn- : 0>4! — '4 4- -4|tOiN- 4 ® — 
j rg, o* itoj O' rg!4 ®'®IO* 0 (0 ® rg, 
i ® I n- !N*! ® I® i to -®- jro oa ;«Ni — ;0 O 

• •, •! •; • • «i • •' • 

O OlOiOlO O O'O ojo o o o 

n , i : i 1 1 1 ■ i •;»».• i i t 


jfH®!~« O to!® — I® m|—l- 
® I® i- 4 |»-«|n* | fO oi® •'O o < 

® NV®I® N-.N- N-!® ® ® < 

® ■ co ! r~ i ® to 4 - m eg — o < 
• «[ «| • • • • | • • • 
O'OiOiO O'o O O O O l 
I • I 1 1 1 1 ] I j I I 1 I 1 1 

. ' I I • 

® rgirg'miN-inj ffl ^ CD i 
® O I ■>*’ I • ' <M . to N- O O < 
| o*;to Ot4|®|® o n -t N i 
|r-’N- T- ! ® |m'4 > 1 - « 1 NH . 
| .1 . • • •'•! 
■OioiO,o;olo o o oioi< 
1 1 i 1 * i ; i •, i i i i ; i ! 

I . 11 . ill 


fsj O' to o fsj,r^|® o ® ro O' — o 

rg 4 m ® ® !® (m ,— © im rg ® m 

co rg 4 in 4 rmm >m m ;<n 4 0 to 

® ® n- ® to 4 |m irg >— 1© O — — 

0 O o o o.oio OiOiO o o o 

1 I I I I I I I I I . i 


N- O' O >#• rg ® rg r- m rg 4- — .rg 
® ® — ® mom rg o m oj O fO 
ON N O NVtO N* ® O m ®|OJ 
® N- r- ® ® tnl 4 mrg rg — < olo 


UJ! I/I ’ *| 

U «J CO NI Zi 
O O UJ —i 

Z * 

n <t a. i 
3 

u o 

x Z 

■a: n uj i 
Ol-U •: 
iuu< V Z; 
to U. —j 


.O K-Z 
— :0 O ! 

NlU- • 
O CO N- OJ Z 1 
uj o H 
co D 


im rg < • 

st'U. > Z 
O ® CL — 

'O • 3 

4 O' to 

® 4 


. || i UJ'lO ♦: 

l ii. JI^MZ 1 

> a O O uj — < 

• UJ Z CL 

» n- n. < a 

3 

a o 
X z 

O co OC t- UJ 

® O O t- o • 

O 11. UJ < v z 

O m co u. — 

• • -< ® 

o — — z o 

O i /> 

II z 

II Ol-Z 

- o o 

o t- UJ — • 

a o co t- m z 

or a uj o — 

uj X to o 


,0 o o o o O O •-« co to ® O rg 
®;O l Kto:iO l/t-0®®®4NN 
rg ® — 4 |® ®f»® 0 ' 0 *ON-cg 
O'® ®N-i® ®4 rorgrg-^OO 
• •••«•• •••• •• 
OO O OlO OOOOOOOO 

I - I I I I I I I I I I I I 


rg— «Om-H4mo®(^OfO® 
mooo-iii'wii'OOO'M' 
4 0'*-* rg — ooo — rg404 
® N- N- ®;® 4 ®rg — OO — — 
• •••••••••••• 

o O O OO OOOOOOOO 
I I I I I I I I I I 


® — ® o — * ® — rg4CT®4rg 
® — O o O ITI r- 4 N- JJ fT O' ® 
i^OH'0'^ , 'l(*i^44®®^ 
O' ® N-N ®®4fOrg— OO 
• •••••••••••• 

ooooooooooooo 

I I I I I I I I t I I I I 


o* m tft ®.® it*®®®®— *N* rg 
® ® N* ® ® 4®rg — OOO— • 
®-*4®®®®®®® f ^O'4 
® ® N- ®® 4frtrsi — OOO — 
• • •« ••••••••• 

OOO ©lO OOOOOOOO 

I I I I I I I I I I 

0 'ir'r^f\jrrir«N 0 'rg» 0 f^N-® 
ON“®N> 0 ®® 0 'rn®f^ ®® 
0'4® — ®4trv®®0'rg®rg 

®®N-N®® 4 mrg — — oo 
ooooooooooooo 

I I I I I I I I I I I I I 


OO (TO 4® ^ KtOOMTCC 

CM4f^t44N-4 4®®®^^ 
04 

®®N-®®4f < >rg— OO^— • 
• ••»••••••••• 

OOOOOOOOOOOOO 
I I I I I I I I I I 

N-®®®OrO®rj — ®<T®N- 
O®^®— •O'co4 — 44rO — 

O® — tr\® 0 '<-*rgofnvr'®m 
a co tc k -o in in 4 m rvj — * o O 
• •••••••••••• 

ooooooooooooo 

I I I I I I I I I I I I I 


lT'( 7 '<Njrg — — rg 4 g 3 — irt tn,® 

rr o N m ® <r\ <d o ® 4 O' rg 4 

C 7 > 4 ® ® N- N 1 ® ®lTitn N o,in 

® ® N ® in 4 m rg — O O 

• •«•••••••••■• 

0 o c o o o,o 000 o o;o 

1 I I I I I I I I I 


g’-jN-H'OriOBNONO® 
co r- inmo®®m — (7'— ‘cnm 
Nrg®Nr-®®®®tnrg< 7'4 
®®N-®®40(M — OOO — 

ooooooooooooo 

I I I t I I I I I I 


rg O' ro OO (M®irg,® m rg ® 1 o 
n ® in 4 n- O' >rg in n- o ® m 4 
n rt ao rg 4 ®jo* — • o ® o 4 0 
N r» ® ® in 4 rt irsj — o uiu 
• • • • • •! •• • *j • • 

0 0 OOO O 1 O 1 O o o o ojo 

1 1 1.1 1 1111:1 1 1,1 


:m O' rg 4 o ®i4 m ® rg ®.4 m 
i® ® O O O 014 OiflO.ai® mlrg 
j® rg 4 ® 4 4jmtm rg]rai4 *-t|® 

!®.® ‘N-i® itn.4 m|rg 1— 1 O 1 O.— — 

i« • I • 

i o o 0 0 O'OiO O 0:0,0 0 0 

I t I I I II tl I , I 


o — N- ® tn‘0 4i0' ® 4 ®iO 

O' ® m O' O' ® o O IN- rg rg O' ] 4 

: ® in — in ® — 4 ® in* O' rg ml— 

N- n N- ® in in 4 rn rg — — 0 0 

'•♦••••••••••• 

000000 0000 o 00 
II I I I till'll I I ill 


O —I 

QL CL 
UJ UI II 
>- N II 

z z 

UI UI O' 
OO CL 

CL O 
UJ UJ UJ X 
-J -J -j o 

o o 

CL CL 


in rg o o olo 000 ® tn(o 

o o — m rn 4i® ® r* ® ® o > iO 

• •••o •:••••••!• 

00000 oto 00 000 — 


125 


LER = 0.0060 TER = 0.0060 LER = 0.0060 TER = 0.0060 LER = 0.0060 TER » 0.0060 

CHORD = 1.3306 CHORD = 1.3435 CHORD = 1.3521 



Table C-l. Stage E Airfoil Coordinates (Continued) 
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CIRCLE CENTER -C.E832 -0.8016 -0.9721 -0.8119 -0.9576 -0.8255 

CIRCLE CENTER -0.0936 0.1571 -0.0538 0.1563 -0.0550 0.1555 

LER = 0.0060 TER = 0.0060 LER = 0.0060 TER = 0.0060 LER - 0.0060 TER * 0.0060 

CHORD = 1.3572 CHORD = 1.3560 CHORD = 1.3552 



Table C-l. Stage E Airfoil Coordinates (Continued) 
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Table C-l. Stage E Airfoil Coordinates (Continued) 
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LE CENTER C. C 223 -0.0577 0.0140 -0.0626 0.0104 -0.0649 

LE CENTER C.6179 0.8759 0.9083 O.R205 0.9176 0.8126 

LEW = 0.0061 TER = 0.0088 LER = 0.0061 TFR = 0.0066 LER = 0.0061 TER = 0.0064 

CHORD = 1.2350 CHORD = 1.2695 CHORD = 1.2747 



Table C-l. Stage E Airfoil Coordinates (Continued) 
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Table C-l. Stage E Airfoil Coordinates (Continued) 
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LER = 0.0061 TER = 0.0115 LER = J.0060 TER = 0.C077 LER - 0.0060 TER = 0.0059 

CHORD = 1.2060 CHORD = 1.2109 CHORD = 1.2168 



Table C-l. Stage E Airfoil Coordinates (Continued) 
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LE CENTER 1.2004 -0.C839 1.2011 -0.0817 1.2024 -0.0683 

LER = 0. 0060 ItR = 0.0060 LER = 0.0060 TER = 0.0060 LER = 0.0060 TER = 0.0060 

CHORD = 1.2085 CHORD = 1.2093 CHORD = 1.2116 



Table C-l. Stage E Airfoil Coordinates (Continued) 
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